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Abstract 

This article has a 3-fold objective: i) to provide a panoramic view of several types 
of 3D vs 4D approaches in Field Theory (Tamm-Dancoff, Bethe Salpeter Equation 
(BSE), Quasi-potentials, Light-Front Dynamics, etc) for strong interaction dunam- 
ics; ii) to focus on the role of the Markov- Yukawa Transversality Principle (MYTP) 
as a novel paradigm for an exact 3D-4D interlinkage between the corresponding 
BSE amplitudes; iii) Stress on a closely parallel treatment of qq and qqq BSE's 
stemming from a common 4-fermion Lagrangian mediated by gluon (vector)-like 
exchange. The two-way interlinkage offered by MYTP between the 3D and 4D BSE 
forms via a Lorentz-covariant 3D support to the BS kernel, gives it a unique status 
which distinguishes it from most other 3D approaches to strong interaction dynam- 
ics, which give at most a one-way connection. Two specific types of MYTP which 
provide 3D support to the BSE kernel, are considered: a) Covariant Instantaneity 
Ansatz (CIA); b) Covariant LF/NP ansatz (Cov.LF). Both lead to formaly identi- 
cal 3D BSE reductions (thus ensuring common spectral predictions), but their 4D 
manifestations differ sharply: Under CIA, the 4D loop integrals suffer from Lorentz 
mismatch of the vertex functions, leading to ill-defined time-like momentum inte- 
grals, but Cov LF is free from this disease. Some practical uses of MYTP as a basis 
for evaluating various types of 4D loop integrals are outlined. 
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1 Introduction: Effective BSE-SDE Framework 



Ever since the success of the Tomonaga-Schwinger-Feynman-Dyson formalism in QED [1], 
corresponding field-theoretic formulations have been in the forefront of strong interaction 
dynamics since the early fifties, the main strategy being to device various 'closed' form 
approaches which are represented as appropriate 'integral equations'. One of the earliest 
efforts in this direction was the Tamm-Dancoff formalism [2] which showed a great intuitive 
appeal. In this method, the state vector of the system under consideration is Fock- 
expanded in terms of a complete set of eigen- functions of the free field Hamiltonian, so 
that the expansion coefficients are the successive "amplitudes for finding present in the 
field a definite number of bare particles with definite spins and momenta...." [3]. This 
method was first systematically applied by Dyson (+ Cornell collaborators) in the early 
fifties, to the meson-nucleon scattering problem, for a dynamical understanding of the 
'Delta' and other low-energy resonances; (see ref. [3]). It leads to 3D integral equations 
connecting amplitudes for successively higher numbers of meson (and/or nucleon-pair) 
quanta, much as the familiar 4D Schwinger- Dyson equations of QED connect (via Ward 
identities) vertex amplitudes of successively higher orders [4]. 

1.1 3D Reduction of BSE: Quasipotentials, Etc 

The 3D Tamm-Dancoff equation (TDE) [3] and the 4D Schwinger-Dyson equation (SDE) 
[4] have been the source of much wisdom underlying the formulation of many approaches 
to strong interaction dynamics. To these one should add the Bethe-Salpeter equation 
(BSE) [5], which is an approximation to SDE for the dynamics of a 4D two-particle 
amplitude, characterized by an effective (gluon-exchange-like) pairwise interaction, on 
the lines of a " Bethe Second Theory" of the Fifties for the effective N-N interaction, but 
now adapted to the quark level. Although not a fundamental (first principle) approach, 
such as the Chiral Schwinger Theory (CST) of (1 + 1)D, [see the Article by P.P.Srivastava 
in this Book], it has attracted more attention in the contemporary literature (as the 4D 
counterpart of the Schroedinger equation) than any other comparable approach. Perhaps 
it is a fair assessment that there is a degree of complementarity between first principle 
(emphasizing theoretical foundation) approaches like CST and second principle (stressing 
applicational aspects) approaches like BSE, in the sense that the lacunae of one are made 
up by the other, so that both together hold the key to a resolution of the strong coupling 
problem. In this Article we are concerned with the latter only. 

A major bottleneck for the BSE approach has been its resistance to a probability 
interpretation, due to its 4D content. This has led to many attempts at its 3D reductions 
[6-9]: Instantaneous approximation [6]; Quasi- potential approaches [7,8]; variants of on- 
shellness of the associated propagators [9]. In [7,9], the starting BSE is 4D in all details, 
including its kernel, but the associated propagators are manipulated in various ways to 
reduce the 4D BSE to a 3D form as a fresh starting point of the dynamics; in [8], the old- 
fashioned 3D perturbation theory is reformulated covariantly to give a 3D quasipotential 
equation. These methods are briefly sketched in Sect 2. 

At the 4D level, the BSE [5] is still the most widely used form of 2-particle dynamics, 
though the problems of probabilistic interpretation were the chief reason for these 3D 
formulations. Nevertheless the regular 4D equations of the full-fledged SDE-BSE types 
have been widely employed [10], as prototypes of strong interaction dynamics, addressing 
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issues of gauge and chiral symmetries, as well as dynamical breaking of chiral symmetry 
via an NJL-type mechanism [11]. In such full-fledged field-theoretic approaches, the NJL- 
mechanism of contact interaction must of course be replaced by space-time extended 
interaction, and the dynamical breaking of chiral symmetry (DBxS) corresponds to the 
use of SDE for the self-energy operator [12]. As a general remark, while for conceptual 
issues impinging on formulational self-consistency, there is little alternative to the full- 
fledged 4D equations, their applications to physical systems must recognize some ground 
realities. For example, the mass spectra of hadrons (which are revealed in Nature as 
0(3)-like [13]), suggest that the role of the time-like dimension (although an integral part 
of the dynamics) is not on the same footing as that of the space-like dimensions, so that 
a naive expection of 0(4)-like spectra [14] may be quite misleading. Indeed this issue is 
quite central to the very theme of this article, viz., 3D-4D interlinkage of BS-amplitudes, 
and will claim attention throughout. 

An alternative form of 2-particle dynamics (which also contributes to reducing the 
effective degrees of freedom from 4D to 3D) is that of Dirac constrained Hamiltonian 
formalism [15], developed by Komar and others [16]. The logic of this approach is that 
constraints Hi have a twin role, viz., they not only 'constrain the motion' in phase space, 
but also generate it in their 'Hamiltonian' capacity. These 'constraints' must be mutually 
compatible in the sense [Hi, Hj] = 0. Such compatibility relations restrict the dependence 
of the interaction on the relative time t, and require a 'reciprocity relation' between the 
constituent potentials, something akin to Newton's III Law. Such descriptions are valid 
for both two-boson and two-fermion systems [16], in the sense of coupled Klein-Gordon 
and Dirac equations respectively. As this formalism is reviewed in detail in this Book by 
Lusanna [17a], it will not be discussed further here. In the same spirit, more fundamental 
approaches like the Chiral Schwinger Theory (CST) in (1 + 1)D, may be found in [17b]. 

1.2 Light Front(LF)/ Null Plane(NP) Dynamics 

A powerful form of 3D dynamics came into prominence after Weinberg's discovery that 
the dynamics of the infinite momentum frame [18] serves as a cure for many ills in the 
theory of current algebras, by greatly simplifying the rules of calculation of Feynman 
diagrams of old fashioned perturbation theory. In the present context of strong inter- 
action dynamics, the great virtue of Weinberg's infinite momentum method [18] lies in 
the simplicity and transparency of the integral equations for multiparticle potential scat- 
tering problems [19]. Indeed, the structure of the 3-momenta (p±,p\\) appearing in this 
formalism is but a paraphrase for the standard null-plane variables first introduced by 
Dirac to project his theme [20] that a relativistically invariant Hamiltonian theory can 
be based on 3 different classes of initial surfaces (space-like, time-like, and null-like). The 
structure of such a Hamiltonian theory is strongly dependent on these respective surface 
forms whose "stability groups" (i.e., those generators of the Poincare group that leave the 
initial surface invariant) , are 6,6,7 respectively, thus giving the 'highest score' (7) to the 
null-plane dynamics (x = x 3 ) whose 'kinematical' generators form a closed algebra, and 
include among others the quantity P + = P + -P3 (which plays the role of the 'mass' term 
77 in the Weinberg notation [18]). On the other hand, the dual generator P_ = P _ P? is 
the 'Hamiltonian' of the theory. 

Leutwyler and Stern [21] gave a covariant 3D formulation of the Dirac theme [20] 
in terms of null-plane variables. A more explicit covariant formmulation in the null- 
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plane language was given by Karmanov [22a] using diagrammatic techniques with on-shell 
propagators and spurions, on the lines of the Kadychevsky approach [8], which has been 
recently reviewed by Carbonell et al [22b] , (referred to as KK). All these methods, including 
the Wilson group's [23], give rise to 3D integral equations for a (strongly interacting) two- 
body system, bearing strong resemblance to the other 3D BSE forms [6-9] above. Again, 
there is no getting back to the original 4D BSE form, the nearest connection being a 
one-way reduction from 4D BSE to 3D on the covariant null-plane [22b]. 

1.3 Markov- Yukawa Transversality: 3D-4D Inter linkage 

Finally we come to a rather novel approach of more recent origin [24-25], based on the 
Markov- Yukawa Transversality Principle (MYTP) [26]. To motivate this approach, it is 
necessary to go back in time to Yukawa's non-local field theory [26b] according to which 
the field variable is a function of both x and p, unlike in local field theory in which the 
field variable is only a (local) function of x. Although unacceptable for an elementary 
particle/field, the non-local field theory is ideally suited to a composite particle, whose 
extended structure effectively provides for a momentum dependence in the direction of the 
total 4- momentum P^. Indeed the Yukawa theory [26b] was in a way the forerunner of a 
later theory of bi-local fields A4(x, y) [27] for the formulation of the Effective Action for a 
2-body dynamical system [27]. This approach was employed by the Prevushin group [25] in 
their formulation of the relativistic Coulomb problem in the Salpeter approximation [6b] in 
a covariant form, with the choice of the preferred direction governed by the 4-momentum 
P M of the composite as the canonical conjugate to its cm. position X = (x + y)/2: 
P = —idx- More specifically the MYTP [26] is expressed by the condition [25]: 

Z ^ M (^X) = 0; z = x-y (1) 

where the direction P M guarantees an irreducible representation of the Poincare' group for 
the bilocal field M [26c]. This condition in turn is equivalent to a covariant 3D support 
to the input 4-quark Lagrangian, whence follow the SDE and BSE as equations of motion 
with 3D support to the effective BSE kernel under covariant instaneity. 

Alternatively, the 3D support ansatz may be directly postulated at the outset for the 
pairwise BSE kernel K [24] by demanding that it be a function of only g M = q — q.PP^/P 2 , 
which implies that q.P = 0. In this approach, the propagators are left untouched in their 
full 4D forms. This is somewhat complementary to the approaches [6-9] (propagators 
manipulated but kernel left untouched), so that the resulting equations [24-25] look rather 
unfamiliar vis-a-vis 3D BSE's [6-9], but it has the advantage of allowing a simultaneous 
use of both 3D and 4D BSE forms via their interlinkage. Indeed what distinguishes the 
Covariant Instantaneity Ansatz (CIA) [24] from the more familiar 3D reductions of the 
BSE [6-9] is its capacity for a 2-way linkage: an exact 3D BSE reduction, and an equally 
exact reconstruction of the original 4D BSE form without extra charge [24]: the former 
to access the observed 0(3)-like spectra [13], and the latter to give transition amplitudes 
as 4D quark loop integrals [24]. (In the approach of the Pervushin Group [25], however, 
the built-in 3D-4D interconnection which follows from MYTP [26], apparently remained 
unnoticed in their final equation). In contrast the more familiar methods of 3D give at 
most a one-way connection, viz., a 4D — > 3D reduction [6-9], but not vice versa. 
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At this point it is perhaps worth noting that even the Salpeter equation [6b] had 
(in principle) the ingredients for a reconstruction of the 4D BS amplitude \1/ in terms of 
3D ingredients, provided the 'instant' form, (see eq.(2.1) in Sect 2), of the interaction 
kernel had been employed on the RHS of the 4D BSE form, and simultaneously the 4D 
BS amplitude \1/ on the RHS had been eliminated in favour of the 3D BS amplitude tp, 
exactly as was done under CIA [24] . This would have amounted to using the Transversality 
Principle [26] (albeit non-covariantly), but this feature had apparently remained unnoticed 
by subsequent workers who continued to employ the Salpeter equation [6b] in its 3D form 
only. 

1.4 Q CD- motivated Effective Lagrangians 

The Transversality Principle (MYTP) [26] underlying the 3D-4D interconnection [24], 
termed 3D-4D-BSE in the following, of course needs supplementing by physical ingredients 
to govern the structure of the BSE kernel, much as a Hamiltonian needs a properly defined 
'potential'. However its canvas is broad enough to accommodate a wide variety of kernels 
which must in turn be governed by independent physical principles. In this respect, short 
of a full-fledged QCD Lagrangian approach, the orthodox view (which we adopt) is to 
stick to an effective 4-fermion Lagrangian as a starting point of the dynamics, from which 
the successive equations of motion (SDE, BSE, etc) follow in the standard manner. [As 
already noted at the outset, this is in keeping with the Bethe Second Principle Theory for 
effective N — N potentials as an input for the physics of the nuclear many-body problem]. 

In particular, a basic proximity to QCD is ensured through a vector- type interaction 
[12], which while maintaining the correct one-gluon-exchange structure in the perturba- 
tive region, may be fine-tuned to give any desired structure to the intermediate gluon 
propagator in the infrared domain as well. Although empirical, it captures a good deal of 
physics in the non-perturbative domain while retaining a broad QCD orientation, albeit 
short of a full-fledged QCD formulation. More importantly, the non-trivial solution of 
the SDE corresponding to this generalized gluon propagator [12] gives rise to a dynamical 
mass function m(p) [12] as a result of DBxS, w.r.t. an input Lagrangian whose chiral 
invariance stems from a vector-type 4-fermion interaction between almost massless u — d 
quarks. These considerations form the standard basis for a Lagrangian-based BSE-SDE 
framework [10] for Dynamical Breaking of Chiral Symmetry (DBxS) [11], for a space- 
time extended 4-quark Lagrangian mediated by vector exchange [12]. This generates a 
mass- function m(p) via Schwinger-Dyson equation (SDE), which accounts for the bulk of 
the constituent mass of ud quarks. The same BSE-SDE formalism [12,10] can be simply 
adapted [28] to the MYTP [26] -based 3D-4D-BSE formalism [24] which reproduces 3D 
spectra of both hadron types [29] under a common parametrization [28] for the gluon 
propagator, with a self-consistent SDE determination [28] of the constituent mass; see 
Sect 3. 

A BSE-SDE formulation [10] on QCD lines represents a 4D field- theoretic general- 
ization of 'potential models' [30], wherein the generalized 4-fermion kernel [12] represents 
the non-perturbative gluon propagator, which can be easily adapted [28] to MYTP [26]). 
The 4D feature of BSE-SDE gives this framework a ready access to high energy ampli- 
tudes, while its 'off-shell' features give it a natural access to hadronic spectra [13]. It has 
thus an interpolating role between (low energy) quarkonia models [30], and (high energy) 
QCD-SR [31] techniques whose domains are largely complementary; details may be found 
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in a recent review [32]. 

1.5 MYTP: Cov Instantaneity vs Cov LF/NP 

While MYTP [26] ensures 3D-4D interconnection [24] under covariant instantaneity ansatz 
(CIA) in the composite's rest frame [24], its main disadvantage lies in the ill-defined na- 
ture of 4D loop integrals which acquire time-like momentum components in the expo- 
nent ial/gaussian factors associated with the different vertex functions, due to a 'Lorentz- 
mismatch' among the rest-frames of the participating hadronic composites. This problem 
is especially serious for triangle loops and above, such as the pion form factor, while 2- 
quark loops [33] just escape this pathology. This problem was not explicitly encountered 
in the light-front (LF/NP) ansatz [34] in an earlier study of 4D triangle loop integrals, 
but this approach was criticized [35] on grounds of non-covariance. The CIA approach 
[24] which made use of MYTP [26], was an attempt to rectify the Lorentz covariance 
defect, but the presence of time-like components in the gaussian factors inside triangle 
loop integrals, e.g., in the pion form factor [36], impeded further progress. 

In an attempt to remedy this situation, a generalization of MYTP [26] was proposed 
recently [37] to ensure formal covariance without having to encounter time-like compo- 
nents in the gaussian wave functions appearing inside the 4D loop integrals. The desired 
generalization was achieved by extending the Transversality Principle [26] from the co- 
variant rest frame of the (hadron) composite [24], to a covariantly defined light-front [37] 
(Cov LF). It was found that while preserving the 3D-4D BSE interconnection, the result- 
ing 3D equation under Cov LF [37] turns out to be formally identical to the old-fashioned 
null-plane formalism [34,38], so that the latter enjoys ipso facto covariance (despite its 
'looks'). This 'covariant' LF/NP method [37] stands fairly direct comparison with other 
covariant LF approaches [22-23]. 

1.6 Scope of the Article 

This article has a 3-fold objective: A) a bird's eye view of some principal 3D vs 4D 
dynamical methods for the strong interaction problem that have been proposed over 
the last half century; B) Putting in perspective a novel property of the Markov- Yukawa 
Transversality Principle (MYTP), viz., a 2- way 3D-4D interconnection in the BS dynamics 
of 2- and 3-quark hadrons; C) Stressing a close parallelism between qq and qqq BSE's which 
stem from a common 4-fermion Lagrangian mediated by gluon (vector)-like exchange. 
Especially noteworthy is the capacity of MYTP [26] to achieve a 3D-3D interlinkage, a 
property which has remained obscured from view in the contemporary literature, vis-a- 
vis more familiar approaches to BSE and allied forms of dynamics [6-10, 18-23] which 
are either 3D or 4D in content, but have no provision for any interlinkage between these 
two dimensions. While the details of individual MYTP applications (several of which 
are dealt with in [32]), are not a part of this Article, the practical uses of MYTP in the 
strong interaction dynamics regime will nevertheless be a focus of attention by virtue of its 
distinct advantage in the evaluation of 4D loop integrals with arbitrary vertex functions, 
while providing easy access to the spectroscopy sector. To that end, an outline of the 
dynamical structure of some principal 3D methods, [7-9], [18-23], is provided in Sect. 2 as 
a background for comparison on 4D loop integral techniques, while on the issue of hadron 
spectra, which are basically 0(3)-like [13], a comparison between non-MYTP [7-9] and 
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MYTP [24-26] forms of dynamics does not bring out new physics. 

For a better understanding of the working of MYTP, it will be necessary to present 
two types in parallel for comparison, viz., Covariant Instantaneity or CIA [24-25], and 
Covariant Light-front (Cov LF) [37], which demand that the BSE kernel K for pairwise 
interaction be a function of relative momentum q transverse to the composite 4-momentum 
in the first case [24], and to the Covariant Null-plane in the second [37]. It will be shown 
that both types lead to identical 3D BSE reductions (so that their spectral predictions are 
formally the same), but their reconstructed 4D vertex functions reveal profound differences 
in structure: The Lorentz mismatch of individual wave functions that characterizes the 
CIA form [24], leading to complex amplitudes [34], disappears in the alternative Cov LF 
approach [37], but in general such integrals are dependent on the light-front orientation n M , 
as in other covariant approaches [22-23] . To eliminate such terms, a simple prescription of 
'Lorentz completion' seems to suffice to produce an explicitly Lorentz invariant quantity 
such as was shown for the pion form factor [37]; (alternative prescriptions exist in other 
LF/NP formulations [22b]). For a historical perspective, it is useful to recall that in the 
old-fashioned NPA approach too [38], a very similar result had been found for various 
types of triangle loop amplitudes [36], despite a lack of manifest covariance [35] in that 
approach, but now MYTP [26] on the covariant null-plane [37] fills this formal gap. 

1.7 Outline of Contents 

Sect. 2 briefly outlines some historical approaches to an effective 3D form of strong inter- 
action dynamics: Levy-Salpeter [6]; Logunov- Tavkhelidze [7a]; Blankenbecler-Sugar [9]; 
Todorov [7d] ; Weinberg [18]; Feynman et al [39]. Sect. 3 provides the theoretical frame- 
work with a short derivation of the BSE-SDE from an input chirally invariant Lagrangian, 
incorporating the original CIA form [25,24] of MYTP [26], on the lines of ref.[25] in terms 
of bilocal fields [27]. It also includes a derivation [28] of the dynamical mass function 
m(p) for an understanding of the constituent mass via Politzer additivity [40]. Sect. 4 
collects some basic results on the null-plane formalism due to Leutwyler- Stern [21], es- 
pecially the role of the 'Angular Condition' in ensuring a formal 0(3)-like invariance. 
From Sect. 5 onwards, the focus is on MYTP [26]-orientation for bringing out its unique 
property which distinguishes it from most other approaches, viz., the 3D-4D interlinkage 
of BS amplitudes. 

Sect. 5 gives a comparative view of the working of MYTP on the BSE forms in CIA 
[24] versus Cov LF [37], and outlines the derivation of the 3D BSE, as well as an explicit 
reconstruction of the 4D BS wave function in terms of 3D ingredients, with 3-momentum 
Q — (q±, 93)5 where the third component emerges as a P-dependent one, suitably adapted 
to the CIA [24] or Cov LF [37] respectively. Sec. 6 gives a corresponding derivation for the 
3D-4D interconnection for a qqq BSE structure under CIA conditions. Sec. 7 illustrates, 
through the calculation of triangle-loop integrals, the relative advantage of Cov LF [37] 
over the CIA [24] version of MYTP [26], in producing a well-defined structure for the 
pion e.m. form factor in a fully gauge invariant manner, and illustrating in the process 
the method of 'Lorentz-completion' for explicit Lorentz invariance, with the expected k~ 2 
behaviour at high k 2 . MYTP also gives a more general structure of triangle loop integrals 
for three-hadron form factors. Sect. 8 summarises our conclusions. 
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2 Quasipotentials And Other 3D Dynamical Equa- 
tions 



The reduction of the 4D BSE for an N — N pair to the 3D level in the Instantaneous 
Approximation was first investigated in the non-adiabatic domain of pseudoscalar meson 
theory (effect of pair-creations included) by Levy [6a], who showed that this 3D BSE 
form is entirely equivalent to the corresponding Tamm-Dancoff equations [2] in the same 
(non-adiabatic) limit. 

On the other hand, Salpeter [6b] employed the adiabatic approximation (no pair cre- 
ation effects) to give a systematic 3D reduction of the fermionic BSE, using projection 
operators for large and small components. The adiabatic approximation amounts to re- 
placing the propagator A F (x — x') for the exchanged meson by 

/inf 
A F (x-x',a;o - x' )d(x - x' Q ) (2.1) 
- inf 

and simply gives the Yukawa potential between two particles. Similarly, in the Instanta- 
neous (adiabatic) Approximation, IA for short, the 4D wave function ty(x) = \l/(x, t) for 
relative motion of two particles becomes simply \l/(x, 0). In the momentum representation, 
these statements read respectively as 

A F (k) A F (k); V(q) = / %)#(q, go) (2.2) 
The Salpeter 3D BSE in the IA for a relativistic hydrogen-atom is [6b] : 

2 

(E - ^(q) - tf 2 (q)) X (q) = / <?k * _ [A 1+ A 2+ - A!_A 2 _] X (k) (2.3) 

where the 3D wave function x(q) is related to the corresponding 4D quantity by an 
equation of the form (2.2), and the symbols A± are energy projection operators for the 
large/small components, etc. 

2.1 Logunov-Tavkhelidze Quasipotentials 

A different form of 3D reduction of the 4D BSE was proposed by Logunov- Tavkhelidze 
[7a] in the language of Green's functions (G-fns) for 2-particle scattering whose momen- 
tum representation my be written as G(pip 2 ; p[p 2 ) (with indicated 4-momenta before and 
after), which satisfies a 4D BSE [7a]: 

(27r) s A(p 1 )A(p 2 )G(p 1 p 2 ;p 1 p 2 ) = 5{pi - p' 1 )5(p 2 - p' 2 )+ J dp'[dp' 2 ^{p 1 p 2 \p'[p' 2 ^G{p'[pl;p l p 2 ) 

(2.4) 

where A(p,) = p\ + mf, etc. Expressing this equation in cm. (P) and relative (g) 
4-momenta, and taking out the <5-fns due to the cm. motion, this equation simplifies 

(2tt) 4 A( Pi ) A(p 2 )G(q, q'; P) = 5(q -?') + / dq"K(q, q")G(q", q; P) (2.5) 

Next, they defined the 3D G-fn for the relative motion as a double integral w.r.t. the two 
time-like momenta: 

G(q,q';P) = J q J ' q' G(q, q'; P) (2.6) 
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Now in operator notation, the 4D BSE (2.5) may be written as G = Go + G$KG, from 
which the kernel K has the formal representation K— Gq 1 — G^ 1 . The L-T trick [7a] now 
consists in using the double integrals on the time-like momenta as in eq.(2.6) to formally 
define the 3D quasipotential K as 

K = Go' 1 - G- 1 (2.7) 

which can be expanded perturbatively in the symbolic form [7a] 

K = Go'GoKGoGo 1 - Go'GqKGoKGoGo 1 - .... (2.8) 

to any desired order of accuracy; [note that the inverse G-fns are just the self-energy 
operators]. If V(q,q';E) is the quasi- potential to a given order of accuracy, then, the 
BSE satisfied by the 3D BS wave function ip(q) is of the form [7a]: 

(E 2 -q 2 - m 2 )^(q) = J d 3 q'V(q, q'; E)^(q') (2.9) 

where the 'denominator function on the LHS arises from integrating Go= A(p 1 ) _1 A(j9 2 ) _1 
w.r.t. q Q and rearranging. 



2.1.1 Narrow resonances in charged particle systems 

Within the last decade, the L-T theory [7a] has witnessed some interesting applications 
[41] to the understanding of 'new' narrow e + e~ resosances observed in heavy ion collisions 
[42]. To that end, the authors [41] have employed an equation of the form (2.9) which 
reads for this system as [41]: 

2c(M - 2ui)m = ^ / ,' 3PW) r (2.10) 

V jrK ^ J (2tt) 3 J 2uj'q{M - u - u' - q + iO) V ; 



where u = ^m 2 + p 2 , and q —\ p — p' |. The results indicate a possible interpretation 
of the observed peaks [42] as new quasi-stationary levels arising from the solution of the 
quasi- potential equation. More interestingly, they also suggest a close relationship of the 
observed states [38] with the von Neumann- Wigner [43] levels embedded in the continuum. 



2.2 Blankenbecler-Sugar Equation 

Another type of quasipotential was proposed by Blankenbecler-Sugar [9], as follows. The 
2-particle scattering amplitude T(q, q') due to a 4D potential V(q, q') in the ladder ap- 
proximation satisfies the BSE [9]: 

T(q,q') = -t(2nr 4 J d 4 q"V(q, g")[m 2 + (P/2+g") 2 ]- 1 [^ 2 + ( J P/2-g' / ) 2 ]- 1 ^(9"^') (2.11) 

where the 2-particle 'free' G-fn is exhibited as the product of the two propagators inside 
the integral on the RHS. To express this equation in 3D form, the B-S [9] trick consists 
first in putting q' on the energy shell, which means that q' Q = 0. and q' 2 = s/4 — m 2 , 
where s = —P 2 is the square of the cm. energy. Next, the on-shell part E2 of the free 
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2-particle G-fn is obtained by taking only the 5-in parts of the two propagators which 
gives rise only to two-particle cuts in the physical region 

E 2 (q") = 2ttJ ds'(s' - sy^lm 2 + (P'/2 + q") 2 ]5[m 2 + (P'/2 - q") 2 } (2.12) 
where s' = —P' 2 , and P' has only a fourth component. This works out as 

E 2 (q") = ^(q'^iq^ + mW 2 - ? 2 )]-l (2-13) 

The balance R 2 of the free G-fn is not singular along the positive cut of the s-variable. If 
it is neglected in the first approximation, and only E 2 from (2.13) is substituted in (2.11), 
then after a trivial integration over q'g, the resultant 3D equation has the form 

T(q q') ~ V(a </) + -/ H ' X^W^ (2 14) 

A comparison of (2.9) and (2.14) shows that although both equations are formally 3D 
in looks, there is a vast difference in their contents: The L-T [7a] form (2.9) involves 
only 3-momenta q = q, since the Hilbert space has been 'truncated' by integrating out 
over their fourth components. The B-S [9] form (2.14) on the other hand has 4-momenta 
formally throughout (no truncation of Hilbert space), except that they are on their mass 
shells ! Thus formal covariance is violated in both equations, although in different ways. 



2.3 Kadychevsky-Todorov Equation 

Still another form of 3D (Lippmann-Schwinger) equation was given by Todorov [7d], 
following the Covariant method of Kadychevsky [8]. In the Todorov approach [7d], the 
potential V w is defined as an infinite power series in the coupling constant which fits the 
perturbative expansion of the scattering amplitude T w for two particles of masses mi,m 2 
and 4-momenta pi,p 2 and q±,q 2 before and after respectively. The quantity T w in the 
off-shell regime satisfies the L-S equation [7d] 

T w {p, q + Mp, q + / d d k— — : 2.15 

ir 2 w J k 2 — b 2 — te 

where the 3D quantities in the cm. frame are defined as 

Pi = -P2 = p; qi = -q2 = q (2 .16) 

and on the energy shell, the corresponding time-like quantities are 

Pio + P20 = w = qi + q 20 ; -p 2 = -q 2 = w 2 ; Aw 2 b 2 (w) = \(w 2 , raj, m\) (2.17) 

This equation too has strong resemblance to the L-T equation [7a]. The corresponding 
equation for the bound state wave function 0(p) is 

7r 2 w{k 2 - b 2 {w))(/){p) = - J d 3 kV{p, k)0(k) (2.18) 

Both B-S [9] and Todorov[7d] equations have been extensively employed in the liter- 
ature. 
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2.4 Infinite-Momentum Frame: Weinberg Equation 



Weinberg [18] observed some remarkable simplifications that occur when the results of 
old-fashioned perturbation theory are expressed in a reference frame in which the total 
3- momentum P is very large. In this limit, the 3- momentum p n of the n-th particle may 
be projected parallel and perpendicular to P, and the results collected as follows: 

p n = r? n P + q n ; q n .P = 0; ^7^ = 1; ^ q n = 0. (2.19) 

n n 

The quantity r\ n > in this theory, plays the role of 'mass' of the n-th particle (in a 3D 
Schroedinger-type equation), and in the P — > oo limit, the rules of calculation become 
very simple: all old-fashioned perturbative diagrams passing through negative energy 
intermediate states vanish, while for the contributing diagrams, the propagator for an 
intermediate state c in a transition from a to b, has the form 2[s a — s c + ie]" 1 , where s for 
any state is the usual total cm. energy squared: 

S = IK + m rJA?n; S a = S b = S c , etc. (2.20) 

n 

Momentum conservation at each vertex is 3D in content: 

^(A^^A^q) (2.21) 

in accordance with the conservation of rj and q, eq.(2.15). The Weinberg counterpart of 
the L-T [7a], B-S [9] and Todorov [7d] equations (2.9), (2.14) and (2.18) respectively, is 
the integral equation [18] 

<q, | T | qr] >=< qr] \ V \ qr] > + J d q J drj _ ^ _ q//2 _ m2 - ~ ] 

(2.22) 

Although this equation is effectively 3D, and has considerable similarity to the correspond- 
ing equations of [7a,7d,9] above, there is a big difference, viz., the angular momentum is 
no longer a well-defined concept in this 3D description. This gap was bridged later by 
Leutwyler-Stern [21] by invoking the 'angular condition' [21], after it became clear that 
the Weinberg approach is equivalent to Dirac's [20] null-plane dynamics; see Sect. 3 below. 



2.5 The FKR Model For 2- And 3-Quark Dynamics 

Before ending this Section, we draw attention for historical reasons, to a unique paper 
by Feynman and collaborators [39], FKR for short, which gave an integrated view of 2- 
and 3-quark hadron dynamics, and played a big role in shaping the direction of strong 
interaction physics to come. The importance of the FKR approach stems, among other 
things, from the fact that these authors were the first to show the way to a unified 
treatment of both 2- and 3-quark hadrons within a common dynamical framework, which 
was to serve as a model for the future. This paper effectively incorporated all the relevant 
aspects of quark dynamics that had been generated piecemeal in the Sixties, and had by 
and large come to be accepted, viz., the group structure 577(6) x 0(3), the symmetrical 
quark model, and harmonic oscillator classification of hadron states (based on their linear 
M 2 -N plots) on the one hand [44], and the mechanism of single-quark transitions, quark 
recoil effects, etc, on the other [45]. 
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The FKR model, which made essential use of harmonic confinement, sought to give a 
relativistic meaning to the internal motion of light quarks through the following definitions 
of 2- and 3-quark Hamiltonians [39, 38]: 

- K M = 2(pl + pi) + -^ 2 { Xl - x 2 f + Const = P 2 + M 2 M ; (2.23) 
-K B = 3( Pl 2 + p 2 2 + p 3 2 ) + 4^ 2 E (*i - x *f + Const = ( p2 + M b) ( 2 - 24 ) 

^ b 123 

where xi^—idi^; p\ = Pi^Pi^, etc. The quantity Q which is postulated to be the same 
for both systems, has the significance of the universal Regge slope lGeV 2 ) as observed 
[44] in their respective spectra; [Note the geometrical factors as cofficients in front of the 
respective kinetic and potential terms above]. The operators K M l B are the 'free' propaga- 
tors (albeit with h.o. confinement) for the mesons and baryons, whose 'poles' correspond 
to the eigenvalues (spectra) of their squared masses. The presence of a perturbation 5K 
can be simulated in a standard gauge- invariant manner, by the substitutions p — > p M — eV^ 
or Pn — > Pn — g^A^ for vector and axial vector couplings respectively, after rewriting p 2 ^ 
as (7-p) 2 , while the i — > j transition amplitudes are just < hj \ 5K \ hi >, by standard 
rules of quantum mechanics. 

A major achievement of the FKR model was its success in giving two distinct types of 
unification, viz., a common framework for Spectroscopy and transition amplitudes; and 
ii) a unified dynamical treatment of qq and qqq hadrons. Both these features represented 
landmarks in a dynamical understanding of the quark model, yet the FKR model failed on 
another count: the 'wrong' sign of the time-like momenta in the gaussian wave functions 
for the hadrons was a disease which pointed to an asymmetric role of time-like (ID) 
momenta vis-a-vis the space-like (3D) ones. Attempts to cure this disease by a Euclidean 
treatment (via Wick rotation) [46a] failed on the spectroscopy front [13] which reveal only 
0(3)-like spectra, while other non-covariant treatments [46b] were not very successful 
either. Nevertheless the lessons from the FKR model were significant pointers to the need 
to treat the ID time-like and 3D space-like d.o.f.'s on different footings in a future quest 
for a covariant theory [24-26]. 

3 Self-Energy And Vertex Fns Under MYTP 

As a first step towards introducing the MYTP [26] theme, we collect in this Section some 
essential machinery for the interconnection between self-energy and vertex functions via 
Schwinger-Dyson (SDE) and Bethe-Salpeter (BSE) equations, starting from a chirally 
invariant Lagrangian characterized by a vector-type interaction [12] as a prototype for 
a gluon-exchange propagator in the non-perturbative QCD regime [28]. To that end, 
we shall first outline the method of bilocal fields [27] to derive the equations of motion 
(SDE and BSE) from such a Lagrangian, following the Pervushin Group's [25] bilocal 
field method, under MYTP [26] conditions of covariant instantaneity [24]. This will be 
followed by a general result connecting self-energy and pion- quark vertex functions in 
the chiral limit, i.e., when the pion mass vanishes. This result in turn paves the way 
to a derivation [28], under MYTP [26] conditions of Covariant Instantaneity [24], of the 
mass function m(p) whose low momentum limit is the main contributor to the constituent 
mass, via Politzer additivity [40]. 
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3.1 Method of Bilocal Fields for BSE-SDE 

The effective action for a system of two interacting massless fermions constrained by 
MYTP [26] is given by 



W e ,M,$\ = J d A x[^x)(i 1 d-m Q )i>(x) + l -J (PyMyMx))K(2+,X)Mx)ij>(y))] (3.1) 

where z = x — y; X = (x + y)/2. z 1 - is the component of z transverse to the P-direction. 
Now redefine the action (3.1) in terms of bilocal fields M via the Legendre transformation 
[27] on the second term to give 

- l - J d 4 xd 4 yM{x, y)K~ 1 {z ± , X)M(x, y) + J d 4 xd 4 y{*p{x)ip{y)M{x, y) (3.2) 
Then in an obvious short-hand notation [25b], (3.1) may be written as 

W eff [M] = (#, (-G, 1 + M)- l -(M, K-'M) (3.3) 

where Go is the inverse Dirac operator for the free fermion field. After quantization over 
N c fermion fields and normal ordering, the action takes the form [25b] 

1 inf 1 

W ef f[M] = -MM, K~ l M) + iN c ]T -<£ n (3.4) 

^ 71=1 71 

where $ = GqM is a matrix in (x, y) space, and its successive powers are defined in the 
standard matrix fashion. Now for the quantization of the action (3.4), its minimum is 
given by 

The corresponding 'classical' (lowest order) solution for the bilocal field is S(x — y) which 
depends only on the difference x — y due to the translation invariance of the vacuum 
solutions. Next expand the action (3.4) around the point of minimum M = S + M', and 
denote the small fluctuations M as a sum over the complete set of 'classical' solutions Y. 
Then in the next order of extremum, we have: 

jj^p \ M '=o r - (3.6) 

Eqs. (3.5-6) give respectively the SDE for E and BSE for T: 

E(x-y) = m 5 4 (x-y) + tK(z ± ,X)Gz(x-y); (3.7) 
T = iK(z ± ,X) J d A Zl d A z 2 Gx(x - z 1 )F(z 1 , z 2 )G^(z 2 - y) 

which describe the spectrum of the fermions and composites respectively. In momentum 
space these equations for the mass operator and vertex function are 

E(p) =m + tj -0^V(p - qhvGdlhm V, = PJ\P\ (3.8) 

m = * J 704^ - ^)7^[G S (g + P/2)Y^)G^q - P/2)] 7? 7 (3.9) 

where Gs(g) = (7g — X^g- 1 )) -1 ; V is the scalar part of the kernel K with 3D support; k is 
the transverse part of k w.r.t. the direction of the total 4-momentum P^. 
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3.2 Self-Energy vs Vertex Fn in Chiral Limit 

The formal equivalence of the mass-gap equation (3.8) and the BSE (3.9) for a pseu- 
doscalar meson in the chiral limit [12] will now be demonstrated by adapting them to a 
non-perturbative gluon exchange propagator [28] with an arbitrary confining form D(k) 
(not just the perturbative form k~ 2 ). The SDE, eq.(3.8), after replacing the color factor 
A1.A2/4 by its Casimir value 4/3, and a relabelling of symbols [28], now reads 

E(p) = ^(2tt)- 4 / d i kD(k)Tf,S' F (p - kht (3.10) 
S' F is the full propagator related to the mass operator E(p) by 

E(p) + i-f.p = S^ip) = A(p 2 )\i 1 .p + m(p 2 )] (3.11) 

thus defining the mass function m(p 2 ) in the chiral limit m c = 0. In the same way, 
eq.(3.9) for the vertex function of a qq hadron (H) of 4-momentum P^ made up of 
quark 4- momenta p 12 = P/2 ± q reads as 

T H (q,P) = -^(2tt)- 4 | d A q'D{q - q'hr}S F (q' + P/2)T H (q', P)S F (q' - P/2) 7? ) (3.12) 

The complete equivalence of (3.10) and (3.12) for the pion case in the chiral limit P M — > is 
easily established. Indeed, with the self-consistent ansatz Th —75^(0), eq.(3.12) simplifies 
to 

T(g) = ^(27r)- 4 | tfhyriS'pik - q)T(q - k)S' F {q - k)^D{k) (3.13) 

where the replacement q' = q — k has been made. Substitution for S' F from (3.11) in 
(3.13) gives 

F(p) = -^(2tt)- 4 / d A k— — ^fe^j — (3.14) 

Kyj 3 V ; J A 2 (p-k){m 2 {{p-k) 2 ) + {p-k) 2 ) v ; 

where we have relabelled q — > p. On the other hand substituting for S' F (3.11) in (3.10) 
gives for the mass term of S(p) the result 

^Wrt = -;i(2-)" 4 / A'^#¥^3T (3-15) 
VF ; VF ; 3 V ; 7 A 2 (q')(m 2 (q' 2 ) + q' 2 ) V ; 

where q' = p — k. A comparison of (3.14) and (3.15) shows their equivalence with the 
identification T(q) = A(q)m(q 2 ) , i.e. the identity of the vertex and mass functions in the 
chiral limit, provided A = 1, (which corresponds to the Landau gauge; see [32]). Although 
obtained here in the context of MYTP [26] this result is independent of this ansatz. A 
more explicit gauge theoretic derivation of the equations for the self-energy and vertex 
functions is given in [32]. 



3.3 Dynamical Mass As DB X S Solution of SDE 

We end this Section with the definition of the 'dynamical' mass function of the quark 
in the chiral limit (M n = 0) of the pion-quark vertex function T(q), in the 3D-4D BSE 
framework [24,28]. The logic of this follows from the BSE-SDE formalism outlined above, 
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eqs.(3. 10-15), for the connection between eq.(3.15) for m(p) and eq.(3.14) for F(q) in the 
limit of zero mass of the pseudoscalar. So, setting M = in the (unnormalized) Hqq 
vertex function T H this quantity may be identified with the mass function m(p), in the 
limit P M = 0, where p M is the 4- momentum of either quark; (note the appearance of the 
'hatted' momentum). The result is [28,32] 

9 aO 
fYl -|- p 

m(p) = [uj(p);\ / 2p.n} q — 0(p) (3.16) 

under CIA and Cov LF respectively. The normalization is such that in the low momentum 
limit, the constituent ud mass m q is recovered under CIA [28], while the corresponding 
'mass' under Cov LF is p + [32] . 

4 Null-Plane Preliminaries 

The Weinberg infinite momentum method received a major boost through an understand- 
ing of Bjorken scaling [47a] in deep inelastic scattering, as well of the Feynman parton 
picture [47b] in the same process. The similarity of the P — > inf and the null-plane de- 
scriptions became clear with the demonstration by Susskind [48] of the U(2) structure of 
the LF/NP language wherein the role of 'mass' is played by the combination p + = Po+p 3 , 
and subsequently a more complete formulation of null-plane dynamics by Kogut and Soper 
[49] within the Hamiltonian formalism in the context of field theory. 

In a different direction, efforts were made to extend the Lorentz contraction ideas 
to finite momentum frames, designed to bring out the effect of Lorentz contraction on 
cluster form factors as a result of motion [50]. In the respect, the role of the Breit frame is 
particularly interesting since it gives the best possible overlap [50b] between the initial and 
final clusters. The Lorentz contraction factors [50] in turn are the key to an understanding 
of 'dimensional scaling' [51], especially in a 'symmetrized' version [50d] of the Breit frame 
[50b] which exactly reproduces the correct 'power counting' [51]. And the Weinberg result 
[18] is duly recovered in the P — > inf limit, giving rise to the null-plane dynamics; (for 
more details of these results, see [38]). A more complete formulation of LF/NP dynamics, 
albeit with a finite number of degrees of freedom was given by Leutwyler- Stern [21] which 
comes closest to the original Dirac [20] spirit, and is summarized below. 

4.1 Leutwyler-Stern Formalism 

Leutwyler-Stern [21] employed a Hamiltonian approach for investigating the properties 
of a relativistic 2-body system with a finite number of d.o.f.'s, and postulating a general 
criterion of 'covariance' in the form of an operator relation among the mass and spin 
operators of the system. Their formalism is based on the maximum 'stability group' (in 
the Dirac [20] sense) for the null-plane surface X0 + X3 = 0, which gives rise to the following 
seven 'kinematical' generators [38]: 

K, = (P 1 ,P 2 ,P + ,E 1 ,E 2 ,K 3 ,J 3 ); 2P ± = P ±P 3 ; (4.1) 
2E 1 = K x + J 2 , 2E 2 = K 2 - J i; Ji = ^e ijk M jk] K t = M 0i . 
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The matrix elements of /C form a closed algebra (r, s = 1, 2): 

[X 3 ,£ r ] = ~iE r ; [K 3 ,P + ] = -iP + ; [J 3 , E r ] = +ie rs E s ; (4.2) 
[J3, P r ] = -\-ie rs P s \ [E r) P^\ = —i5 rs P + 

On the other hand, there are 3 'Hamiltonians' which can be chosen in one of several 
different ways. To that end, it is necessary to introduce certain rotation operators l{ 
defined by U \ n >= Ji \ n >, on a rest system | n >, but extended to states | p, n > 
defined by 

I p, n >= Expi-ifcE! - t(3 2 E 2 - t(3 3 K 3 ] \ n >; f3 r = p r /p + ; f3 3 = \n(2 P+ /M) (4.3) 
such that I commutes with the algebra of /C. More explicitly, 

Ii\p,n> = Expl-i/3^ - i(3 2 E 2 - i(3 3 K 3 ]Ji | n >; (4.4) 

In particular, 

I 3 = J 3 + (E 1 P 2 - E 2 P 1 )/P + = (W + W 3 )/P + ; W, = l -e, vaP P"M^. (4.5) 

is the Pauli-Lubanski operator, and MI r = W r — P r W + /P + , where r = 1,2. One 
thus arrives at the 'dynamical group' D = (M, Jj), or (M 2 , M/j), which has the structure 
of U(2) [48], because of (3.4). D is really a 3-member group, since ^3 already belongs 
to /C by virtue of (3.5). For a particular significance of I, note its connection with the 
non-relativisic Galilei-invariant algebra generated by the momentum P and Galilei boosts 
K, viz., 

I = J + m~ 1 KxP; (m = mass) (4.6) 

Now Galilei invariance of a system is equivalent to the condition that the corresponding 
dynamical algebras constitute a unitary representation of U(2). In the relativistic case, 
there is a superficial similarity to the U(2) structure of D, but unlike the NR case, only 
the component I 3 is now 'kinematical', by virtue of (3.5), while ii i2 are 'dynamical', and 
do not have explicit representations by themselves. L-S [21] sought to bridge this gap 
by imposing the 'covariance' requirement in the form of an 'angular condition' for the 
operators Jj as follows: 

xiMh + x 2 MI 2 + x L I 3 ; x L = P x x x + p 2 x 2 + P + x_ (4.7) 

which can be shown to be valid on the null-plane x + = xq + x 3 = [38] . 

The L-S formalism [21] provides a compact support to the longitudinal momentum z 
of 2-particle system with constituent 4-momenta pi j2 : 

2zP + =p 1+ -p 2+ - P + =p 1+ +p 2+ ; ~<x<+ 1 -. (4.8) 

The internal wave function is defined by 

< Pl , p 2 I P, >= 2P+5 3 ( Pl + p 2 - P)0(q±, x) (4.9) 
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where is a matrix in helicity space (h', h"), with the norm: 



< I >= \ I x £ | Mcu, z) | 2 (4.10) 



4 Z h'/i" 



The L-S [21] structure formally allows the introduction of a 3-vector q and angular mo- 
mentum L for the internal motion of a composite of mass M with (equal) constituent 



masses m q as 



with the phase space 



q=(qi,Mi); L = -iq X V, (4.11) 



= JK» = 4K + q>). (412) 

With these definitions of q and L, the theory formally preserves the concept of L- 
invariance of a qq system despite the apparently asymmetric treatment meted out to 
the transverse and longitudinal components of 3-momenta in the NP or P = inf formal- 
ism [18, 25, 48-9]. This invariance can be traced to angular condition (4.7). Incidentally, 
the historic FKR model [39], despite its other defects, was found by L-S [21] to satisfy 
the angular condition (4.7). 

An alternative but more pedigogical recipe to achieve the same end was given in [38] 
via the simpler condition P. q = 0, to be consistently imposed between the total (P M ) and 
relative (g M ) 4-momenta of a qq system, as outlined in subsection 4.2 below [38]. 



4.2 An Alternative "Angular Condition" P.q = 

For unequal masses of the (quark) constituents with 4-momenta p± t 2, the total (P) 
and relative (q) 4-momenta are given by the Wight- Gaer ding definitions [52] 

2 2 

Pi,2 = m 1 , 2 P±q; 2m 1>2 = l± l — 2 2 (4.13) 

where M = \J — P 2 is the composite mass. The condition Pq = is satisfied on the mass 
shells m\ 2 + p\ 2 = of the respective constituents, by virtue of the Wight man- Gaer ding 
definition (4.13). 

To link the condition P.q = with the construction of an effective 3-vector in the 
null-plane language, so as to preserve the invariance of the angular momentum concept, 
note that this condition translates to the relation q_ = —q + M 2 /P 2 which expresses the 
component g_ in terms of q + , in a frame P± = 0, since in this frame, P+P- = M 2 
on the mass shell of the composite. [The collinearity condition is not a restriction for 
a two-body system]. This relation then allows a definition of the 3-momentum q with 
the components (q_i_>93)> with ^3 = Mq + /P + , which preserves the meaning of L in the 
sense of L-S [21], together with NP covariance. For any other internal 4- vector for the 
composite, a similar 3-vector A may be defined as (Aj_, A3), with A 3 = MA + /P + , via the 
condition AP = 0. Examples of A^ are polarization vectors, Dirac matrices, etc. Using 
these techniques, null-plane wave functions of the L-S type [21] have been constructed and 
applied to hadronic processes via quark loops [34]. A more formal mathematical basis 
for this prescription comes from MYTP [26] on the covariant null-plane [37]; see Sect. 5 
below. 
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5 3D-4D BSE Under MYTP: Scalar /Fermion Quarks 



We now come to our objectives (B) and (C), viz., 3D-4D interlinkage of BS amplitudes 
brought about by MYTP [26], and a unified treatment of qq [24,37] and qqq [53] systems 
under MYTP conditions. The full calculational details with 4-fermion couplings via glu- 
onic propagators have been collected in a recent review [32]. However, to bring out the 
basic mathematical structures, we shall derive the 3D-4D interconnection with spinless 
quarks for 2- and 3-quark systems in this and the next sections respectively. Further, 
we shall consider two MYTP methods in parallel for comparison: i) Covariant Instan- 
taneity Ansatz (CIA) [24-25]; ii) Covariant LF/NP Ansatz (Cov LF) [37], to bring out 
the structural identity of the resulting BSE's for a qq system. This will be followed by 
a reconstruction of the 4D BS vertex functions for both types [24, 37] as basic building 
blocks for 4D quark loop calculations. 

5.1 3D-4D BSE Under CIA: Spinless Quarks 

For a self-contained presentation, with unequal mass kinematics [24], let the quark con- 
stituents of masses mi^ and 4-momenta pi )2 interact to produce a composite hadron of 
mass M and 4-momentum P M . The relative 4-momentum g M is related to these by 

Pi,2 = rh 1>2 P ± q\ P 2 = -M 2 ; 2m lt2 = l±{m\-ml)/M 2 (5.1) 

These Wightman-Garding definitions [52] of the fractional momenta mi i2 ensure that 
q.P = on the mass shells mf+pf = of the constituents, though not off-shell. Now define 
Qn — 9/z — Q-PP/i/P 2 as the relative momentum transverse to the hadron 4-momentum P M 
which automatically gives q.P = 0, for all values of q^. If the BSE kernel K for the 2 
quarks is a function of only these transverse relative momenta, viz. K = K(q,q'), this is 
called the "Cov. Inst. Ansatz (CIA)" [24] which accords with MYTP [26]. For two scalar 
quarks with inverse propagators A 12 , this ansatz gives rise to the following BSE for the 
wave fn $(g,P) [24]: 

i(2n) 4 A 1 A 2 ^(q, P)= J d 4 q'K(q, q')$(q', P); A 1>2 = m\ 2 + p\ 2 (5.2) 

The quantities m 1}2 are the 'constituent' masses which are strictly momentum dependent 
since they contain the mass function m{p) [12,28], but may be regarded as constant for 
low energy phenomena: m(p) = m(0). Further, under CIA, m(p) = m(p), a momentum- 
dependence which is governed by the DB\S condition [28] (see below). 

To make a 3D reduction of eq.(5.2), define the 3D wave function 0(g) in terms of the 
longitudinal momentum Ma as 

(j)(q) = J Mda$(q,P); Ma = Mq.P/P 2 (5.3) 

using which, eq.(5.2) may be recast as 

i(27r) 4 A 1 A 2 $(g,P) = J d z q'K(q,q')<P{q'); d 4 q' = d 3 q' Mda' (5.4) 

Next, divide out by A X A 2 in (5.4) and use once again (5.3) to reduce the 4D BSE form 
(5.4) to the 3D form 

(2vfD{qm) = J <W?,$W); ^ = f%£- 2 ( 5 - 5 ) 
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Here D(q) is the 3D denominator function associated with the like wave function 0(g). 
The integration over da is carried out by noting pole positions of A 12 in the a-plane, 
where 

Ai, 2 = ^ lj2 2 -M 2 (m lj2 ±a) 2 ; w 1)2 2 = m 2 2 + g 2 (5.6) 
The pole positions are given for A lj2 = respectively by 

M(<7 + mi) = ±ui =F ie; M(cr - rh 2 ) = ±u 2 =F *e (5.7) 

where the (±) indices refer to the lower/upper halves of the a- plane. The final result for 
D(q) is expressible symmetrically [24]: 

D{q) = M u D Q (q)\ W = ~ + - ^ 

= f _ M^M5. A = M 4 -2M 2 K + m 2 ) + K-m 2 ) 2 (5.9) 

The crucial thing for MYTP[26] is now to observe the equality of the RHS of eqs (5.4) and 
(5.5), thus leading to an exact interconnection between the 3D and 4D BS wave functions 
[24]: 

r(j)^A 1 A 2 $( 5 ,P) = ^M« (5.10) 

Eq.(5.10) determines the hadron-quark vertex function F(q) as a product D(f> of the 3D 
denominator and wave functions, satisfying a relativistic 3D Schroedinger-like equation 
(5.5). The simultaneous appearance of the 3D form (5.5) and the 4D form (5.4), leading 
to their interconnection (5.10), reveals a two-tier character: The 3D form (5.5) gives 
the basis for making contact with the 3D spectra [13], while the reconstructed 4D wave 
(vertex) function (5.10) in terms of 3D ingredients D and can be used for 4D quark-loop 
integrals in the standard Feynman fashion. Note that the vertex function T = D(f)/(2iir) 
has a general structure, independent of the details of the input kernel K. Further, the 
-D-function, eq.(5.8), is universal and well-defined off the mass shell of either quark. The 
3D wave function is admittedly model- dependent, but together with D(q), it controls 
the 3D spectra via (5.5), thus offering a direct experimental check on its structure. Both 
functions depend on the single 3D Lorentz-covariant quantity q 2 whose most important 
property is its positive definiteness for time-Ike hadron momenta (M 2 > 0). 



5.2 Cov LF/NP for 3D-4D BSE: Fermion Quarks 

As a preliminary to defining a 3D support to the BS kernel on the light-front (LF/NP), 
on the lines of CIA [24], a covariant LF/NP orientation [37] may be represented by the 
4- vector n M , as well as its dual n M , obeying the normalizations n 2 = n 2 = and n.n = 1. 
In the standard NP scheme (in euclidean notation), these quantities are n = (001; —i) / \[2 
and n = (001; i)/y/2, while the two other perpendicular directions are collectively denoted 
by the subscript _L on the concerned momenta. We shall try to maintain the n-dependence 
of various momenta to ensure explicit covariance; and to keep track of the usual NP 
notation p± = p ±p 3 , our covariant notation is normalized to the latter as p + = n.p\^2; 
p_ = —n.p^/2, while the perpendicular components continue to be denoted by p± in both 
notations. 
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In the same notation as for CIA [24] , the 4th component of the relative momentum q = 
m 2 pi — rhiP2, that should be eliminated for obtaining a 3D equation, is now proportional 
to q n = n.q, as the NP analogue [37] of P.qP/P 2 in CIA [24], where P = p 1 + p 2 is the 
total 4- momentum of the hadron. However the quantity q — q n n is still only q±, since 
its square is q 2 — 2n.qn.q, as befits q\ (readily checked against the 'special' NP frame). 
We still need a third component p 3 , for which the correct definition turns out to be [37] 
q 3fM = zP n n^, where P n =P.n and z = q.n/P.n, which checks with q 2 = q\ + z 2 M 2 . We 
now collect the following definitions/results: 

q± = q — q n n; q — q± + xP n n; x — q.n/P.n; P 2 = — M 2 ; (5.11) 
q n ,P n = h.(q,P);q.n = q.n; q.n = 0; P±.q± = 0; 

P.q = P n q.n + P.nq n ; P.q = P n q.n; q 2 = q\ + M 2 x 2 

Now in analogy to CIA, the reduced 3D BSE (wave-fn 0) may be derived from the 4D 
BSE (5.2) for spinless quarks (wave-fn $) when its kernel K is decreed to be independent 
of the component q n , i.e., K = K(q,q'), with q = (q±,P n n), in accordance with MYTP 
[26] condition imposed on the null-plane (NP), so that d A q = d 2 q^dq^dq n . Now define a 
3D wave-fn 4>(q) = J dq n Q(q), as the CNPA counterpart of the CIA definition (5.3), and 
use this result on the RHS of (5.2) to give 

*(27r) 4 $(g) = Ar 1 A 2 - 1 J d 3 q'K(q, q')<j>(q') (5.12) 

which is formally the same as eq.(5.4) for CIA above. Now integrate both sides of eq.(5.12) 
w.r.t. dq n to give a 3D BSE in the variable q: 

{2nfD n {q)m = J d 2 q 1 _'dq,'K{q,q') < P{q') (5.13) 

which again corresponds to the CIA eq.(5.5), except that the function D n (q) is now defined 
by 

J dq n A 1 - 1 A 2 - 1 = 2mD~ 1 (q) (5.14) 

and may be obtained by standard NP techniques [38] (Chaps 5-7) as follows. In the q n 
plane, the poles of A 12 lie on opposite sides of the real axis, so that only one pole will 
contribute at a time. Taking the A 2 -pole, which gives 

2qn = _^ ^±Sb^L (5 . 15 ) 

m 2 P.n — q.n 

the residue of Ai works out, after a routine simplification, to just 2P.q = 2P.nq n + 2P n q.n, 
after using the collinearity condition P±.q± = from (5.11). And when the value (5.15) 
of q n is substituted in (5.14), one obtains (with P n P.n = —M 2 /2): 

D n (q)=2P.n(q 2 - X{M2 4 ™l m2) ); q 2 = ql + M 2 x 2 ; x = q.n/P.n (5.16) 
Now a comparison of (5.12) with (5.13) relates the 4D and 3D wave-fns: 

2m$(q) = D n (q)A 1 - 1 A 2 - 1 cf ) (q) (5.17) 
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as the Cov LF counterpart of (5.10) which is valid near the bound state pole. The BS 
vertex function now becomes T = D n x <j>/(2ni). This result, though dependent on the 
LF/NP orientation, is nevertheless formally covariant, and closely corresponds to the 
pedagogical result of the old LF/NP formulation [38], with D n D + . 

A 3D equation similar to the covariant eq.(5.13) above, also obtains in alternative 
LF formulations such as in Kadychevsky-Karmanov [22b] (see their eq.(3.48)). However 
the independent 4- vector n M (which has no counterpart in [22b]), makes this a manifestly 
covariant 4D formulation without need for explicit Lorentz transformations [22b]. The 
'angular condition' [21] is also trivially satisfied by the effective 3-vector q^ appearing in 
the 3D BSE (5.13). A more important contrast from other null-plane approaches is that 
the inverse process of reconstruction of the 4D hadron-quark vertex, eq.(5.17)), has no 
counterpart in them [22-23], as these are basically 3D oriented, not vice versa. 

For fermion quarks with gluonic propagators, the MYTP formulation needs no new 
principles, except for certain technical details involving slight modifications [54] of the 
BSE structure for easier handling; see [32] for detailed steps. The full 4D wave function 
\I/(P, q) may be expressed as a 4x4 matrix [38,32]: 

= S F ( Pl )T(q) lD S F (-p 2 ); T(q) = N H [l;P n /M\D(q)<j>(q)/2iir (5.18) 

where 7^ is a Dirac matrix which equals 75 for a P-meson, 27^ for a V-meson, 27^75 
for an A-meson, etc. The factors in square brackets stand for CIA and Cov LF values 
respectively. Nh represents the hadron normalization. 

6 The qqq BSE: 3D-4D Interlinkage 

We now come to the aspect of MYTP [26] that governs the inter-relation of 3D and 4D 
Bethe-Salpeter amplitudes for 3-body (qqq)-sy stems, in keeping with a perceived 'duality' 
between meson (qq) and baryon (qqq) systems which necessitates a parallel treatment 
between them. In this respect a fairly comprehensive review of baryon dynamics as a 
3-body relativistic system with full permutation symmetries in all relevant degrees of 
freedom [55] has been given recently [32]. These include: A detailed correspondence 
[56] between qqq and quark-diquark wave functions; Complex HO techniques for the qqq 
problem [57]; fermionic qqq BSE with the same gluon propagator for pair qq interactions 
[29] as employed for qq systems [28] , except for reduction by half due to the color factor; 
and Green's function methods for 3D reduction of the 4D BSE form, plus reconstruction of 
the 4D qqq wave function [53], on the lines of the qq problem [24]. Within the formalistic 
scope of this Article however, we shall merely dwell on the last item, viz., Green's fn 
techniques [53] for a 3D reduction of the 4D BSE, plus reconstruction of the 4D wave 
function, for a qqq system for three identical spinless quarks, keeping in forefront the 
issue of connectedness [58] in a 3-particle amplitude whose signal is the absence of any 
S -function in its structure; (for a detailed perspective, see [32]). 

6.1 Two-Quark Green's Function Under CIA 

As a warm up to the method of Green's functions (G-fns), we first derive the 3D-4D inter- 
connection for the corresponding G-fns for 2-particle scattering of two identical spinless 
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particles, before moving on to the 3-body problem in the next 2 subsections. For sim- 
plicity we shall consider the G-fns near the bound state poles, so that the inhomogeneous 
terms may be dropped. In the notation and phase convention of Section 5, the 4D qq 
Green's fn G(pip 2 ;pi'p 2 ) near a bound state satisfies a 4D BSE (no inhomogeneous term): 

i(27r) 4 G( PlP2 ; Pl W) = A 1 ~ 1 A 2 ~ 1 J dp 1 "dp 2 "K(p 1 p 2 ;p 1 "p 2 ' / )G( Pl 'W;PiP2'); (6.1.1) 
where 

A 1 = Pl 2 + m q 2 , (6.1.2) 

and m q is the mass of each quark. Now using the relative 4- momentum q = (pi — P2}/2 
and total 4-momentum P = pi + p 2 (similarly for the other sets), and removing a 8- 
function for overall 4-momentum conservation, from each of the G- and K- functions, 
eq.(6.1.1) reduces to the simpler form 

i(27c) A G(q.q') = Af'Af 1 J dq"Mda"K(q, q")G(q", q') (6.1.3) 

where % = % — °~P^-, with a = (q.P)/P 2 , is effectively 3D in content (being orthogonal 
to Pfj). Here we have incorporated the ansatz of a 3D support for the kernel K (inde- 
pendent of a and a'), and broken up the 4D measure dq" arising from (6.1.1) into the 
product dq"Mda" of a 3D and a ID measure respectively. We have also suppressed the 
4-momentum P^ label, with (P 2 = — M 2 ), in the notation for G(q.q'). 
Now define the fully 3D Green's function G(q, q') as [53] 

G{q,q') = J J M 2 dada' G(q,q') (6.1.4) 
and two (hybrid) 3D-4D Green's functions G(q,q'), G(q,q') as 

G(q, q') = J MdaG(q, q'); G(q, q') = J Mda'G(q, q'); (6.1.5) 
Next, use (6.1.5) in (6.1.3) to give 

*(27r) 4 G(g, q>) = A 1 ~ 1 A 2 - 1 J dq"K(q, q")G(q", q') (6.1.6) 
Now integrate both sides of (6.1.3) w.r.t. Mda and use the result 

J MdaA l - 1 A 2 ~ 1 = 2mD- 1 (q); D(q) = 4oo(lj 2 - M 2 /4); C? = m 2 + q 2 (6.1.7) 

to give a 3D BSE w.r.t. the variable q, while keeping the other variable q' in a 4D form: 

(2nfG(q, q') = IT 1 / dq"K(q, q")G(q", q') (6.1.8) 

A comparison of (6.1.3) with (6.1.8) gives the desired connection between the full 4D 
G-function and the hybrid G(q, g')-function: 

2mG(q, q') = D(q) Af 1 A^G^q, q') (6.1.9) 
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Again, the symmetry of the left hand side of (6.1.9) w.r.t. q and q' allows rewriting the 
right hand side with the roles of q and q' interchanged. This gives the dual form 

2niG(q, q') = D(q') A^A^Gfa, q') (6.1.10) 

which on integrating both sides w.r.t. Mda gives 

2niG(q, q') = D^A^A^G^, q)- (6.1.11) 

Substitution of (6.1.11) in (6.1.9) then gives the symmetrical form 

(2m) 2 G(q, q') = D(q) A^A^Giq, q')D{q') A^A^ 1 (6.1.12) 

Finally, integrating both sides of (6.1.8) w.r.t. Mda', we obtain a fully reduced 3D BSE 
for the 3D Green's function: 

(27r) 3 G(g, q') = D-\q j dq" K{q,q")G{q" ,q') (6.1.13) 

Eq.(6.1.12) which is valid near the bound state pole, expresses the desired connection 
between the 3D and 4D forms of the Green's functions; and eq(6.1.13) is the determining 
equation for the 3D form. A spectral analysis can now be made for either of the 3D or 
4D Green's functions in the standard manner, viz., 

G(q, q ') = J2<s>n(q;P)K(q';P)/(P 2 + M 2 ) (6.1.14) 

n 

where <£> is the 4D BS wave function. A similar expansion holds for the 3D G-function 
G in terms of (f) n (q)- Substituting these expansions in (6.1.12), one immediately sees the 
connection between the 3D and 4D wave functions in the form: 

2iri$(q, P) = A 1 - 1 A 2 - 1 D(q)<P(q) (6.1.15) 

whence the BS vertex function becomes T = D x <f>/(2ni) as found in [24]. We shall make 
free use of these results, taken as qq subsystems, for our study of the qqq G-functions in 
subsects.6.2-3. 

6.2 3D BSE Reduction for qqq G-fn 

As in the two-body case, and in an obvious notation for various 4-momenta (without 
the Greek suffixes), we consider the most general Green's function G{piP2p 3 ',P\P2Pz) for 
3-quark scattering near the bound state pole (for simplicity) which allows us to drop the 
various inhomogeneous terms from the beginning. Again we take out an overall delta 
function 5{p\ + P2 + P3 — P) from the G-function and work with two internal 4-momenta 
for each of the initial and final states defined as follows [54b] : 

^3^3 = pi - pi ; 3r] 3 = -2p 3 +pi+p 2 (6.2.1) 

P = Pl+P2+P3=Pl' +P2 +P3 (6.2.2) 

and two other sets £i, rji and £ 2 , V2 defined by cyclic permutations from (6.2.1). Further, as 
we shall consider pairwise kernels with 3D support, we define the effectively 3D momenta 
Pi, as well as the three (cyclic) sets of internal momenta £i,f)i, (i = 1,2,3) by [54b]: 

Pi=Pi- ViP ; ii = & - SiP ; fji - UP (6.2.3) 
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Vi = (P. Pi )/P 2 ; = (P.&/P 2 ; U = {P.r]i)/P 2 (6.2.4) 

v^3s 3 = vi — v 2 ; 3t 3 = — 2u 3 + vi + v 2 (+cyclicpermutations) (6.2.5) 
The space-like momenta pi and the time-like ones v { satisfy [54b] 

Pi + P2 + P3 = ; z/i + z/ 2 + z/ 3 = 1 (6.2.6) 



Strictly speaking, in the spirit of covariant instantaneity, we should have taken the relative 
3D momenta £,,f) to be in the instantaneous frames of the concerned pairs, i.e., w.r.t. the 
rest frames of Pij = pi + Pj] however the difference between the rest frames of P and Py 
is small and calculable [54b], while the use of a common 3-body rest frame (P = 0) lends 
considerable simplicity and elegance to the formalism. 

We may now use the foregoing considerations to write down the BSE for the 6-point 
Green's function in terms of relative momenta, on closely parallel lines to the 2-body case. 
To that end note that the 2-body relative momenta are q^ = (pi—pj)/2 = \/3£ k /2, where 
(ijk) are cyclic permutations of (123). Then for the reduced qqq Green's function, when 
the last interaction was in the (ij) pair, we may use the notation G(£kVk] CkVk), together 
with 'hat' notations on these 4-momenta when the corresponding time-like components 
are integrated out. Further, since the pair £k,Vk is permutation invariant as a whole, 
we may choose to drop the index notation from the complete G-function to emphasize 
this symmetry as and when needed. The G-function for the qqq system satisfies, in the 
neighbourhood of the bound state pole, the following (homogeneous) 4D BSE for pairwise 
qq kernels with 3D support: 

*(27r) 4 G(^; £'?/) = £ A 1 _1 A 2 _1 / dq'( 2 Mda 12 "K(q 12 , V; 6V) (6-2.7) 

123 J 

where we have employed a mixed notation (q i2 versus £3) to stress the two-body nature 
of the interaction with one spectator at a time, in a normalization directly comparable 
with eq.(6.1.3) for the corresponding two-body problem. Note also the connections 

012 = VSs 3 /2; q 12 = y/3£ 3 /2; fj 3 = -p 3 , etc (6.2.8) 

The next task is to reduce the 4D BSE (6.2.7) to a fully 3D form through a sequence 
of integrations w.r.t. the time-like momenta Sj,tj applied to the different terms on the 
right hand side, provided both variables are simultaneously permuted. We now define the 
following fully 3D as well as mixed (hybrid) 3D-4D G-functions according as one or more 
of the time-like £,1] variables are integrated out: 

G(£m£'v)= IJ! j dsdtds'dt'G^-W) (6.2.9) 

which is S'3-symmetric. 

G 3v mi'f}') =11 dt 3 dt 3 'G^r]^W); (6.2.10) 

GatitmN) = 11 ds 3 ds 3 'G(^;ev'); (6.2.11) 

The last two equations are however not symmetric w.r.t. the permutation group S3, since 
both the variables £, r\ are not simultaneously transformed; this fact has been indicated in 
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eqs. (8.2.10-11) by the suffix "3" on the corresponding (hybrid) G-functions, to emphasize 
that the 'asymmetry' is w.r.t. the index "3". We shall term such quantities "S^-indexed" , 
to distinguish them from ^-symmetric quantities as in eq. (6.2.9). The full 3D BSE for the 
G-function is obtained by integrating out both sides of (6.2.7) w.r.t. the st-pair variables 
dsids/dtidt/ (giving rise to an S^-symmetric quantity), and using (6.2.9) together with 
(6.2.8) as follows: 

(27T) 3 G(^') = Y.D-\q 12 ) [ d$l 2 K{q 12 ,$l 2 )G{iY-A'f}') (6.2.12) 

123 J 

This integral equation for G which is the 3-body counterpart of (6.1.13) for a qq system in 
the neighbourhood of the bound state pole, is the desired 3D BSE for the qqq system in a 
fully connected form, i.e., free from delta functions. Now using a spectral decomposition 
for G 

G(imi'v') = J2Miv;Pmi'v';P)/(P 2 + M 2 ) (6.2.13) 

n 

on both sides of (6.2.12) and equating the residues near a given pole P 2 = —M 2 , gives the 
desired equation for the 3D wave function for the bound state in the connected form: 

(2nfm-P) = J2D- 1 (q 12 ) [ dq f ; 2 K(q 12 ,q'; 2 )cl>(i"f,";P) (6.2.14) 

123 

Now the ^-symmetry of in the (£j,?ft) pair is a very useful result for both the solution 
of (6.2.14) and for the reconstruction of the 4D BS wave function in terms of the 3D wave 
function (6.2.14), as is done in the subsect.6.3 below. 

6.3 Reconstruction of 4D qqq Wave Function 

We now attempt to re-express the 4D G-function given by (6.2.7) in terms of the 3D G- 
function given by (6.2.12), as the qqq counterpart of the qq results (6.1.12-13). To that end 
we adapt the result (6.1.12) to the hybrid Green's function of the (12) subsystem given 
by G 3v , eq.(6.2.10), in which the 3-momenta 173,173 play a parametric role reflecting the 
spectator status of quark #3, while the active roles are played by q± 2 , q\ 2 = v^fe, £3')/% 
for which the analysis of subsect.6.1 applies directly. This gives 

(2ni) 2 G 3r] ^ 3 f, 3 ; = D(q 12 ) A^A,' 1 G(i 3 f, 3 ; ^M^V^A^ 1 (6.3.1) 

where on the right hand side, the 'hatted' G-function has full S'3-symmetry, although (for 
purposes of book-keeping) we have not shown this fact explicitly by deleting the suffix 
'3' from its arguments. A second relation of this kind may be obtained from (6.2.7) by 
noting that the 3 terms on its right hand side may be expressed in terms of the hybrid G 3 £ 
functions vide their definitions (6.2.11), together with the 2-body interconnection between 
(£3) £3') an d (£3, £3) expressed once again via (6.3.1), but without the 'hats' on i] 3 and 773'. 
This gives 

(VSm) 2 G(^ m ;^ m ') = (V3m) 2 G(Z m ev') 

= ^Ar 1 A 2 - 1 (« v / 3) / 'dq^Mda^'Kiq^q'^G^V-^rj^ 

123 J 

= Yl D(q 12 )A 1 - 1 A 2 - 1 G 3 ^ 3 r ]3 ; i^A^Ar 1 (6.3.2) 

123 
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where the second form exploits the symmetry between £, r\ and rf . 

At this stage, unlike the 2-body case, the reconstruction of the 4D Green's function 
is not yet complete for the 3-body case, as eq. (6.3.2) clearly shows. This is due to the 
truncation of Hilbert space implied in the ansatz of 3D support to the pairwise BSE kernel 
K which, while facilitating a 4D to 3D BSE reduction without extra charge, does not have 
the complete information to permit the reverse transition (3D to 4D) without additional 
assumptions. To fill up this gap in this "inverse" mathematical problem, we look for a 
suitable ansatz for G 3 g on the RHS of (6.3.2) in terms of known quantities, so that the 
reconstructed 4D G-function satisfies the 3D equation (6.2.12) exactly, as a check-point. 
We therefore seek a structure of the form 

Gszdsm; iW) = G(i 3 m; £&) x F(p 3 ,p 3 ') (6.3.3) 

where the unknown function F must involve only the momentum of the spectator quark 
#3. A part of the rj 3 ,r] 3 ' dependence has been absorbed in the G function on the right, 
so as to satisfy the requirements of S^-symmetry for this 3D quantity [53]. 

As to the remaining factor F, it is necessary to choose its form in a careful manner so 
as to conform to the conservation of 4-momentum for the free propagation of the spectator 
between two neighbouring vertices, consistently with the symmetry between p 3 and p 3 . 
A possible choice consistent with these conditions is: 

F(ps, Ps) = C 3 A 3 - 1 5(z/ 3 - W) (6.3.4) 

Here A 3 _1 represents the "free" propagation of quark #3 between successive vertices, 
while C 3 represents some residual effects which may at most depend on the 3-momentum 
p 3 , but must satisfy the main constraint that the 3D BSE, (6.2.12), be explicitly satisfied. 

To check the self-consistency of the ansatz (6.3.4), integrate both sides of (6.3.2) w.r.t. 
ds 3 ds 3 'dt 3 dt 3 to recover the 3D S3- invariant G-function on the left hand side. Next, in 
the first form on the right hand side, integrate w.r.t. ds 3 ds 3 on the G-function which 
alone involves these variables. This yields the quantity G 3 £. At this stage, employ the 
ansatz (6.3.4) to integrate over dt 3 dt 3 . Consistency with the 3D BSE, eq.(6.2.12), now 
demands 

C 3 j j rfz/ 3 rfz/3 / A 3 - 1 5(z/3 - ia,') = 1; (sincedt = dv) (6.3.5) 

The ID integration w.r.t. dv 3 may be evaluated as a contour integral over the propagator 
A -1 , which gives the pole at v 3 = cu 3 /M, (see below for its definition). Evaluating the 
residue then gives 

C 3 = in/(Mu> 3 ); cu 2 3 = m q 2 + p\ (6.3.6) 

which will reproduce the 3D BSE, eq.(6.2.12), exactly] Substitution of (6.3.4) in the 
second form of (6.3.2) finally gives the desired 3-body generalization of (6.1.12) in the 
form 

3G(^;£V) = E^(?i2)AiFA 2F D(^ 12 )A 1F 'A 2F 'G(| 3 ^;&3)[A3F/(M7r ( I;3)] (6.3.7) 

123 

where for each index, A F = — iA^ 1 is the Feynman propagator. To find the effect of 
the ansatz (6.3.4) on the 4D BS wave function $(£77; P), we do a spectral reduction like 
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(6.2.13) for the 4D Green's function G on the LHS of (6.3.2). Equating the residues on 
both sides gives the desired 4D-3D connection between $ and 0: 

P) = g D( j 12 )Ar' A 2 " V«.i; J>) x f" 3 '^ (6-3.8) 
defines the 4D wave fn in terms of piecewise vertex fns Vi, as 

From (6.3.8-9), we infer the baryon-ggg vertex function V 3 corresponding to the 'last' 
interaction in the 12-pair as 

V 3 = D(q 12 <i>(irj) x yj2A 3 6{u 3 *M* - c2$ (6.3.10) 

and so on cyclically. (The argument of the 5-function inside the radical for V 3 simplifies 
to Ps + rriq). This expression had been obtained earlier from intuitive considerations [54b]. 

To account for the appearance of the ID 5-fn under radical in (6.3.10), it is explained 
elsewhere [53] that it has nothing to do with connectedness [58] as such, but merely reflects 
a 'dimensional mismatch' due to the 3D nature of the pairwise kernel K [24] imbedded 
in a 4D Hilbert space. (For a physical explanation, see [53]). A further self-consistency 
check on (6.3.10), is found by taking the limit of a point interaction, which amounts to 
setting K = Constant, when the radical (expectedly) disappears, and gives a Lorentz- 
invariant result [53], in agreement with the so-called NJL-Faddeev (contact) model [59] 
for 3-particle scattering. For the fermion qqq case with pairwise gluonic interactions, the 
details may be found in [60], wherein the strength of the 'color' qq interaction [29] is half 
of that of qq [28]. For brevity, we skip the MYTP [26] derivation of the 4D qqq vertex 
function under Cov LF [37] conditions, which parallels that for the 2-body case [Sect. 5], 
except for the remark that the old-fashioned LF/NP treatment [38] gives the same results 
as the more formal Cov LF treatment in Sect. 5, so that a similar Cov LF form for qqq 
dynamics should be expected [38], with D(q) — > D n (q), etc in (6.3.10). 

7 Triangle Loops Under MYTP On Cov LF/NP 

In this Section, we shall illustrate the MYTP techniques on the covariant light-front to 
bring out the main feature, viz., structure of the triangle loop integrals free from the 
anomalies of time-like momenta in the product of gaussian vertex functions, such as 
complexities in the pion form factor [36] (see Sects. 1 and 5). To that end, we shall 
mainly consider the mathematical structure of the P-meson form factor, followed by a 
brief stetch of the structure of 3-hadron form factors, in the next few subsections, leaving 
routine calculational details to [37,32]. 
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7.1 Pion Form Factor by Cov LF/NP Method 




Figure 1: Triangle loop for e.m. vertex 



Using fig. 1 above, and an identical one with 1 — > 2, (c.f., figs, la, lb of [34b]), the 
Feynman amplitude for the h — > hi + 7 transition is given by [34b] 

2P,F(k 2 ) = A{2nfN n {P)N n {P')em l J + {1 2]; {7A) 

AT (i) = Tr [ 75 ( mi _ i^.p^i^rri! - 27-^1)75(^2 + n-P2)\; Aj = m- + (7.2) 

Pi, 2 = m h2 P ± q; p[ 2 = m ly2 P' ± q' p 2 = p' 2 ; P - P' = Pl - p[ = k; 2P = P + P'. 

(7.3) 

After evaluating the traces and simplifying, T M becomes 

= (P2» ~ P,)[Sm 2 — M 2 — A 2 ] - k 2 p 2 ,/2 + (A x - A;)A; M /4 (7.4) 

The last term in (7.4) is non-gauge invariant, but it does not survive the integration in 
(7.1), since the coefficient of k^, viz., Ai — A[ is antisymmetric in p 1 and p[, while the 
rest of the integrand in (7.1) is symmetric in these two variables. Next, to bring out the 
proportionality of the integral (7.1) to P M , it is necessary to resolve p 2 into the mutually 
perpendicular components p 2 ±, (p 2 .k/k 2 )k and (p 2 .P/ P 2 )P, of which the first two will 
again not survive the integration, the first due to the angular integration, and the second 
due to the antisymmetry of k — p\ —p[ in p 1 and p[, just as in the last term of (7.4). The 
third term is explicitly proportional to P M , and is of course gauge invariant since P.k = 0. 
(This fact had been anticipated while writing the LHS of (7.4)). Now with the help of 
the results 

p 2 .P = -m 2 M 2 -Ai/4-A' 1 /4; 2m 2 = 1 - {m\ - m\)/M 2 - P 2 = -M 2 - k 2 /4, (7.5) 

it is a simple matter to integrate (7.1), on the lines of Sec. 5, noting that terms proportional 
to Ai A 2 and A' 1 A 2 will give zero, while the non- vanishing terms will get contributions only 
from the residues of the A 2 -pole, eq.(5.15). Before collecting the various pieces, note that 
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the 3D gaussian wave functions 0, 0', as well as the 3D denominator functions D n , D' n , do 
not depend on the time-like components p 2n , so that no further pole contributions accrue 
from these sources. (It is this problem of time-like components of the internal 4-momenta 
inside the gaussian 0-functions under the CIA approach [24], that had plagued a earlier 
CIA study of triangle diagrams [36]). To proceed further, it is now convenient to define 
the quantity q.n = p 2 .n — m 2 P.n to simplify the 0- and D n - functions. To that end define 
the symbols: 

(q,q') = q±m 2 k/2; z 2 = q.n/P.n; k = k.n/P.n; {9 k , rj k ) = 1 ± fc 2 /4 ( 7 - 6 ) 
and note the following results of pole integration w.r.t. p 2n [38]: 

J dp 2n i-[l/A i; l/Ai; 1/(A 1 A' 1 )] = [1/D n ; 1/D' n , 2p 2 .n/(D n D' n )} (7.7) 

Details of further calculation of the form factor are given in [37]. An essential result is 
the normalizer N n (P) of the hadron, obtained by setting k^ — 0, and demanding that 
F(0) = I. The reduced (Lorentz-invariant) normalizer N H = N n (P)P.n/M is given by 
[32,37]: 

Nh 2 = 2M(2tt) 3 J (Pqe-flf [(1 + 5m 2 /M 2 )(q 2 - A/4M 2 ) + 2m x m 2 (M 2 - 5m 2 )} (7.8) 

where the internal momentum q = (q±,Mz 2 ) is formally a 3-vector, in conformity with 
the 'angular condition' [21]. The corresponding expression for the form factor is [32, 37]: 

F{k 2 ) = 2MN 2 H (2itfexp[-(Mm 2 k/ (3) 2 /44](vr/3 2 ) 3/2 ^m 1 G(fc) + [1 2] (7.9) 
where G(k) is a function of k ; see eqs.(A. 12-13) of [32]. 
7.2 'Lorentz Completion' for F(k 2 ) 

The expression (7.9) for F(k 2 ) still depends on the null-plane orientation n M via the 
dimensionless quantity k = k.n/P.n which while having simple Lorentz transformation 
properties, is nevertheless not Lorentz invariant by itself. To make it explicitly Lorentz 
invariant, we shall employ a simple method of 'Lorentz completion' which is merely an 
extension of the 'collinearity trick' empolyed at the quark level, viz., P±.q± = 0; see 
eq.(5.11). Note that this collinearity ansatz has already become reduntant at the level of 
the Normalizer N^, eq.(7.8), which owes its Lorentz invariance to the integrating out of the 
null-plane dependent quantity z 2 in (7.8). This is of course because Nfj depends only on 
one 4-momentum (that of a single hadron), so that the collinearity assumption is exactly 
valid. However the form factor F(k 2 ) depends on two independent 4-momenta P, P', for 
which the collinearity assumption is non-trivial, since the existence of the perpendicular 
components cannot be wished away! Actually the quark- level assumption P±.q± = 
has, so to say, got transferred, via the g-integration in eq.(7.9), to the hadron level, as 
evidenced from the /c-dependence of F(k 2 ); therefore an obvious logical inference is to 
suppose this /c-dependence to be the result of the collinearity ansatz Pj_.P^ = at the 
hadron level. Now, under the collinearity condition, one has 

P.P' = P±.P' ± + P.nP'.h + P'.nP.h = P.nP' n + P'.nP n ; P.h = P n . (7.10) 
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Therefore 'Lorentz completion' (the opposite of the collinearity ansatz) merely amounts 
to reversing the direction of the above equation by supplying the (zero term) P±.P' ± to 
a 3-scalar product to render it a 4-scalar! Indeed the process is quite unique for 3-point 
functions such as the form factor under study although for more involved cases (e.g., 
4-point functions), further assumptions may be needed. 

In the present case, the prescription of Lorentz completion is relatively simple, being 
already contained in eq.(7.10). Thus since P,P'= P±k/2, a simple application of (7.10) 
gives 

k.nk n = +k 2 ; P.nP n = —M 2 — k 2 /4; (7.11) 
4k 2 

k = 4M^ = 4 *'- 4 = 4 - 4 * 

This simple prescription for k automatically ensures the 4D (Lorentz) invariance of F(k 2 ) 
at the hadron level. (For comparison with alternative methods [22b], see [37]). 



7.3 QED Gauge Corrections to F(k 2 ) 

While the 'kinematic' gauge invariance of F(k 2 ) has already been ensured in Sec. 7.1 above, 
there are additional contributions to the triangle loops - figs. la and lb of [34b] - obtained 
by inserting the photon lines at each of the two vertex blobs instead of on the quark lines 
themselves. These terms arise from the demands of QED gauge invariance, as pointed 
out by Kisslinger and Li [61] in the context of two-point functions, and are simulated by 
inserting exponential phase integrals with the e.m. currents. However, this method (which 
works ideally for point interactions) is not amenable to extended (momentum-dependent) 
vertex functions, and an alternative strategy is needed, as described below. 

The way to an effective QED gauge invariance lies in the simple-minded substitution 
Pi — eiA(xi) for each 4-momentum pi (in a mixed p, x representation) occurring in the 
structure of the vertex function. This amounts to replacing each q^ occurring in T(q) = 
D(q)<f>(q), by q^ — e q A^, where e q = rh^ei — mi^, and keeping only first order terms in 

after due expansion. Now the first order correction to q 2 is — e q q.A — e q A.q, which 
simplifies on substitution from eq.(7.11) to 

- 2e q q.A = -2e g A M [g M - q.nh^ + P.nq.nn^/ P.n] (7.12) 

The net result is a first order correction to T(q) of amount e q j(q).A where 

j(q), = -4M>g M 0(g)(l - (q 2 - A/4M 2 )/2/3 2 ) (7.13) 

The contribution to the P-meson form factor from this hadron-quark-photon vertex (4- 
point) now gives the QED gauge correction to the triangle loops, in the form of a similar 
function F\{k 2 ) which works out as [37]: 

F,{k 2 ) = A{2^N 2 H e q m l M 2 > J d'q{M 2 - 5m*W[^j^ + ^£^H + {1^2} 

' . ' . ' ' (7 ' 14) 
where M > = sup{M;m 1 + m 2 } [37], and the common factor 2P^ has been extracted as 

for F{k 2 ) in (7.1). Note that e q is antisymmetric in '1' and '2', signifying a change of sign 
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when {1 =>- 2} is added on the RHS. The pole integration of Fi(k 2 ) now yields a result 
like (7.9) for F(k 2 ); see [37] for details. 

The large and small k 2 limits of F(k 2 ) and Fi(k 2 ) are on expected lines, and we 
summarise only the final results for completeness [37]. For large k 2 , the functions F(k 2 ) 
and Fi(k 2 ) both yield the correct asymptotic form C/k 2 , where C = 0.35GeV 2 , to be 
compared with the experimental value [62a] 0.50 ± 0.10, and the (perturbative) QCD 
value [63] 8na s f 2 = 0.296. 

For low k 2 , on the other hand, an expansion of F, F\ in powers of k 2 yields a value of 
the charge radius R according to < R 2 > = — Vfc 2 (F(/c 2 ) + Fi{k 2 )) in the k 2 = limit. 
Of the two functions, only F(k 2 ) contributes in this limit [37]. The numerical values for 
the kaon and pion radii, vis-a-vis experiment [62b], are 

R K = 0.63fm vs(0.53fm); R n = 0.661/m vs(0.656fm). (7.15) 
7.4 Three-Hadron Couplings Via Triangle-Loops 



F(P) 




Figure 2: 3-hadron coupling 

For a large class of hadronic processes like H — > H'+H" and H — > H'+^y, the quark tri- 
angle loop [64] represents the lowest order "tree" diagram for their evaluation. Criss-cross 
gluonic exchanges inside the triangle-loop are not important for this description in which 
the hadron-quark vertices, as well as the quark propagators are both non-perturbative, and 
thus take up a lion's share of non- perturbative effects. This is somewhat similar to the 
'dynamical perturbation theory' of Pagels-Stokar [65], in which criss-cross diagrams are 
neglected. 

We now indicate in the barest outline, the structure of the 3-hadron loop integral 
for the most general case of unequal mass kinematics mi =^ m 2 ^ m 3 , while referring for 
notational details to ref.[32]. The full structure of the 3-hadron amplitude may be written 
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down from fig. 2 above (c.f., fig. 1 of [64]), just as for the e.m. form factor (7.1): 

exhibiting cyclic symmetry, where the normalized vertex function Tj in CNPA [41] is given 
in an obvious notation by eq.(5.18) as 

m) = N^y^Dmm); a = - mh ^; ml) ) (7.17) 

where the 'reduced' denominator function D; L = D i+ Mi/ ' P i+ and the (invariant) normalizer 
N iH is iVj. The color factor and the effect of reversing the loop direction are given by 
2/v^3, etc [38,64]. the overall BS normalizer [38]. 

To evaluate (7.16), we first write the cyclically invariant measure: 

d% = d p ± ^d(x 2 l )M?dy i ; x i =p i+ /P i+ ; y i =p i -/P i . (7.18) 

The cyclic invariance of [7.18] ensures that it is enough to take any index, say 2, and first 
do the pole integration w.r.t. the y 2 variable which has a pole at y 2 = $, 2 =u} 2± / (M^ x 2 ) ■ 
The process can be repeated, by turn, over all the indices and the results added. Note 
that the ^-functions do not include the time-like yi variables under CNPA [37], so that the 
residues from the poles arise from only the propagators. The crucial thing to note is that 
the denominator functions D\ and D 3 sitting at the opposite ends of the p 2 -\ine (c.f. Fig.l 
of [64]) will cancel out the residues from the complementary (inverse) propagators A3 and 
Ai respectively. Indeed by substituting the pole value y 2 = £2, i n Ai^, the corresponding 
residues in an obvious notation work out as [32]: 

Ai ;2 = 6*132 Ml + x 2 n 23 M 2 - 2/} 21 M 3 2 ; A 3;2 = -&rai 2 A/f - x 2 n 21 M* - 2/i 23 M 2 (7.19) 

It is then found, with a short calculation [32], that 

^3(93) ,, Di(qi) , 7 ,„v 
— = 2M 3 x 2 n 23 ; — = 2M l x 2 n 2l (7.20) 

^l;2 ^3;2 

which shows the precise cancellation mechanism between the Dj-functions and the residues 
of the propagators Aj at the A 2 pole. This mechanism thus eliminates [24, 64] the 
(overlapping) Landau-Cutkowsky poles that would otherwise have caused free propagation 
of quarks in the loops. The same procedure is then repeated cyclically for the other 
two terms arising from the A 31 poles. Collecting the factors, the result of all the 3 
contributions is compactly expressible as [64, 32]: 

A(3H) = 8^S 123 j J M^n^dx^xllTR^D^U^M.N^ (7.21) 

where the limits of integration for both variables are — inf < (£2,2:2) < +inf, since these 
are governed, not by the on-shell dynamics of standard LF methods [22-23], but by off- 
shell 3D-4D BSE. The difference from [64] (under CIA [24]) arises from using CNPA [37] 
which has ensured that the (gaussian) functions <pi on the RHS of (7.21) are now free from 
time-like momenta (unlike in CIA [24,64]). 
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Eq.(7.21) is the central result of this exercise. Its general nature stems from the use 
of unequal mass kinematics at both the quark and hadron levels, which greatly enhances 
its applicability to a wide class of problems which involve 3-hadron couplings, either 
as complete processes by themselves (such as in decay processes) or as parts of bigger 
diagrams in which 3-hadron couplings serve as basic building blocks. What makes the 
formula particularly useful for general applications is its explicit Lorentz invariance which 
has been achieved through the simple method of 'Lorentz Completion' on the lines of 
sect. 7.2 for the e.m. form factor of P-mesons; for more details, see [32]. 

As regards two- quark loops, such as for £77(2) mass splittings of P-mesons [33b], and 
the mixing of p and uj off-shell propagators [33a], the distinction between CIA [24] and 
CNPA [37] is less sharp, (no time-like momentum problems in the overlap integrals). The 
same holds for one-quark loops, e.g., in the problem of vacuum condensates. For a review 
of these processes, as well as for other references, see [32]. 

8 Retrospect And Conclusions 

In keeping with our objectives (A)-(C) defined at the outset, Sects. 1-2 have attempted 
a panoramic view of several standard approaches to 3D BSE reductions [6-9] under the 
general Bethe Second Principle Philosophy of effective quark-pair interaction. In par- 
ticular, the relative unfamiliarity with MYTP [26] in the literature, especially its novel 
feature of effecting an exact 3D reductions of the qq and qqq BSE's, as well as exact 
reconstructions of the 4D amplitudes in closely parallel fashions, have necessitated the 
introduction of some background techniques under one roof. To that end Sect. 3 collects 
a general derivation of the equations of motion in interlinked BSE-SDE form from an 
input 4-fermion Lagrangian for 'current' quarks, under MYTP conditions [25], much like 
the derivation of similar equations [10-12] without this constraint. And in preparation 
for the derivation of MYTP-governed equations in a Covariant LF/NP framework, Sect. 4 
collects some essential background material, especially the 'angular condition' [21], under 
Cov LF/NP conditions. With this background, Sect. 5 outlines a comparative derivation 
of the MYTP-controlled 3D-4D interlinkage of qq Bethe-Salpeter amplitudes under both 
CIA [24] and Cov LF [37] conditions. And in keeping with the parallelism between the 
2- and 3-quark treatments, Sect. 6 gives a similar derivation for the qqq sustem. Now this 
twin facility which does not seem to exist in the other 3D approaches [6-9,22-23], gives 
rise to a natural two-tier description [38] , the 3D BSE form being appropriate for making 
contact with the hadron spectra [13], while the reconstructed 4D BSE yields a vertex 
function which allows the use of standard Feynman diagrams for 4D loop integrals. To 
appreciate why two distinct forms of MYTP [26] have been developed on parallel lines: 
Covariant Instantaneity Ansatz [24] (CIA), and Covariant Light-front [37] (Cov LF/NP), 
the advantage of the latter over the former in producing well-defined triangle loop inte- 
grals has been demonstrated in Sect. 7 through the examples of pion form factor and more 
general 3-hadron couplings, except for a (less serious) problem of dependence on the 'null- 
plane orientation' which can be handled through a simple device of 'Lorentz completion' 
and yields an explicitly Lorentz-invariant structure. Similarly the baryon-quark vertex 
function (6.3.10) is a key ingredient of MYTP, for the calculation of baryonic loop inte- 
grals, of which the baryon self-energy [60] is the simplest example, but the calculational 
details [60] have been omitted for brivity. 
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In keeping with its mathematical (formalistic) emphasis of this Article, we have re- 
frained from discussing the phenomenological applications, but it has been shown that the 
canvas of MYTP [26] is broad enough to accommodate additional physical principles. In 
particular, the physical basis chosen for detailed presentation, has been a QCD motivated 
4-fermion Lagrangian (with an effective gluonic propagator) which generates the BSE- 
SDE structure by breaking its chiral symmetry dynamically (DB\S) [11-12], formulated 
within an MYTP [26] framework. 

Clearly, the MYTP is a very powerful Principle which helps organize a whole spectrum 
of phenomena under a single umbrella. For its applications, only a few examples have 
been indicated, but its potential warrants many more. More importantly, the interlinked 
3D-4D structure of BS dynamics under MYTP [26] premises, gives it access to a whole 
range of physical phenomena, from spectroscopy to diverse types of loop integrals. The 
emphasis on the spectroscopy sector as an integral part of quark physics was first given 
by Feynman et al [39]. 
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